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We show that the same physical mechanism is fundamental for two seemingly 
different phenomena such as the formation of two-level systems in glasses and the 
Boson peak in the reduced density of low-frequency vibrational states g{uj)/uj'^ . This 
mechanism is the vibrational instability of weakly interacting harmonic modes. Below 
some frequency ujc <C ojq (where ujq is of the order of Debye frequency) the instability, 
controlled by the anharmonicity, creates a new stable universal spectrum of harmonic 
vibrations with a Boson peak feature as well as double-well potentials with a wide 
distribution of barrier heights. Both are determined by the strength of the interaction 
I (X LOc between the oscillators. Our theory predicts in a natural way a small value for 
the important dimensionless parameter C = P^"^ j pv^ ~ 10~^ for two-level systems 
in glasses. We show that C {Wjlvjo^"^ oc I~^ and decreases with increasing of 
the interaction strength I. The energy W is an important characteristic energy in 
glasses and is of the order of a few Kelvin. This formula relates the TLS's parameter 
C with the width of the vibration instability region cjc which is typically larger or of 
the order of the Boson peak frequency w;,. Since IiWc > IiWh ^ W the typical value 
of C and therefore the number of active two-level systems is very small, less than 
one per ten million of oscillators, in a good agreement with experiment. Within the 
unified approach developed in the present paper the density of the tunneling states 
and the density of vibrational states at the Boson peak frequency are interrelated. 

PACS numbers: 61.43.Fs, 63.50+x, 78.30.Ly 



I. INTRODUCTION 

One of the most typical low-temperature properties of glasses is the existence of two- 
level systems (TLS's) which determine at low temperatures, typically below a few Kelvin, 
and frequencies less than 1 GHz, phenomena such as specific heat, thermal conductivity, 
ultrasonic and microwave absorption and many others (such as the echo, etc)^. Another 
remarkable universal property of almost all glasses is the Boson peak in the low frequency 
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inelastic scattering, proportional to the reduced density of the vibrational states g{uj) juJ^. 
Compared to the TLS's it is observed at much higher frequencies, between 0.5 and 2 THz, 
and persists to higher temperatures, sometimes up to the glass transition temperature 
Tg. Usually these two important glassy features are considered separately and no definite 
connection between them has been established so far. 

The purpose of the present paper is to propose a physical picture in which these two 
seemingly different phenomena are closely interrelated. We will show that the formation 
of the Boson peak in the reduced vibrational density of states (DOS) inevitably leads 
to the creation of two-level systems (and vice versa). A mechanism of the Boson peak 
formation, implying also the formation of TLS's was proposed in our recent papers^i^. 
This mechanism was based on the phenomenon of vibrational instability of weakly in- 
teracting harmonic oscillators (HO's). The instability takes place in the low frequency 
region, < u; < ujc, below some characteristic frequency Uc which is proportional to 
the interaction strength /. Changing the interaction strength one can vary the width 
of the instability region, the position of the Boson peak and the number of TLS's. The 
weakness of the interaction / implies that Uc ujq, where Uq is of the order of Debye 
frequency. Thus the vibrational instability occurs far below the Debye frequency uo and 
has little influence on the high frequency vibrations. 

The low-frequency harmonic oscillators we are speaking about are realized in glasses 
as quasi-local vihisitions (QLV) which can be understood as /oca/ low frequency vibrations 
bilinearly coupled to the sound waves^"^. The existence of quasi-local vibrations in 
glasses has been confirmed in a number of papers (see e.g. the literature cited in Ref.-). 
The local low frequency vibrations are the cores of the QLV's. The cores represent 
collective low frequency vibrations of small groups of atoms^i^. If one plots the potential 
energy against the amplitude of one of these modes^*^, one obtains a mode potential, 
as treated in the soft potential model^^'-^. The vibrational instability results from the 
weak interaction I of these soft oscillators with high-frequency oscillators (with ~ tuo). 
As a result of this harmonic instability and the anharmonicity of the glass the initial 
vibrational density of states go{i^) is reconstructed in the low frequency region < < 
ujc where instability takes place and shows the universal Boson peak feature. 

The microscopic origin of these high frequency oscillators, whose frequencies are much 
higher than the Boson peak one, is not important in this picture. It has been shown 
that, in general, the frequency of the sound waves can pass through the loffe-Regel limit 
near the Boson peak frequencySii2,iiJ,iiiii5,. Therefore, the higher frequency modes might 
have a rather complex structure. They interact with the soft oscillators. This interaction 
causes the vibrational instability and hence the Boson peak and TLS's formation. 

The vibrational instability is a rather general phenomenon and occurs in any system 
of bilinearly coupled harmonic oscillators. It can be considered in a purely harmonic 
approximation. For example, a system of two oscillators with masses Mi 2 and frequen- 
cies uji^2 becomes unstable if the interaction / between the oscillators exceeds a critical 
valued Ic = uJ\UJ2^ M\M2. If one of two frequencies is small then the critical interaction 
Ic is also small. In our case we have such an instability due to the interaction of low and 
high frequency oscillators. Therefore, we can simplify the picture to treat it analytically. 
The physical reason for the instability in this case is the fact that the high frequency 
modes adiabatically follow the motion of the low-frequency oscillator (a description of 
this fact can be given within the adiabatic approximation) . The squared frequency of the 
soft oscillator is reduced by an amount proportional to the strength of the interaction 
squared and, therefore, even can turn negative. For the example of two interacting os- 
cillators with frequencies lvi -C u)2 the new frequency uji of the low-frequency oscillator 
is given by^ 

^i=coni-iyi!). (1.1) 
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It is important that it decreases to zero linearly with Ic — I {J^l oc Ic — I, when I ^ 1^). 

This is the essence of the vibrational instabihty phenomenon in the general case. As 
a result, switching on the interaction between the oscillators, we shift all low frequency 
modes with < towards the boundary point uj = and some of them will cross 
this point (so that the corresponding uj"^ turns negative). Therefore, in such a case (if 
there is no hard gap around zero in the initial density of states, go{^^)) we will have a 
constant distribution of renormalized u"^ around a; = 0. From that we immediately get 
the right wing of the Boson peak. Indeed, the local anharmonicity does not change this 
important property of uniform distribution of renormalized u"^ around zero. One can 
show^ that it restores the mechanical stability of the system by simply reflecting all the 
negative u"^ values back to the positive u;^ range (like in a mirror but with additional 
stretching factor of 2, which is obviously not important, see Eq. fl2.16p ). The strength 
of the anharmonicity itself plays no role in this mirror transformation. Now a constant 
distribution of u"^ around zero (on the positive side) obviously leads to a universal linear 
law for the density of states, g{uj) oc in the interval < u < Uc (independent of the 
initial DOS gQ{uj)). And in its turn this linear u dependence of the reconstructed DOS 
just gives us the right wing of the Boson peak, since g{(jj)/uj'^ (xl/uj and this dependence 
is also universal and independent of the initial DOS gf){uj). 

We want to stress that this DOS transformation (due to the phenomenon of vibra- 
tional instability) is rather general and universal since any monotonous "traffic" of uJ^ 
from positive to negative values (due to interaction between the oscillators, or due to 
changing the temperature or pressure, etc.) always gives a constant distribution of u"^ 
around zero. Therefore, the universal linear DOS, g{oj) oc u;, and the corresponding 
universal right wing of the Boson peak, g{ijj)/uj'^ oc 1/a;, inevitably emerge as a result of 
this instability. 

If the origin of the right wing of the Boson peak looks somehow natural, the left wing 
of the Boson peak appears as a result of the less obvious additional transformation of the 
linear DOS at smaller frequencies. The point is that when the anharmonicity restores 
the mechanical stability of the system, single-well potentials describing the unstable 
soft modes with negative uj"^ are replaced by double-well ones. This means that the 
effective potential energy of the glass in direction of the local soft mode has two minima 
separated by a rather low barrier. Thus, in this scheme, the two-level systems are 
created simultaneously with the Boson peak due to the same mechanism of vibrational 
instability. Besides their own high importance (the TLS's physics in glasses and beyond) 
these double-well potentials play also an important role in building the left wing of the 
Boson peak. It can be explained as follows. 

Due to bilinear interaction between the oscillators, lijXiXj, double- well potentials, 
with a particle vibrating in one of the wells and therefore having a non-zero average 
displacement Xio ^ 0, create random static forces /j ~ -^y-a^jo acting on other oscillators. 
In a purely harmonic case, these linear forces would not affect the frequencies (and 
linear density of states g{uj) oc uj would not change). However, together with the local 
anharmonicity the static forces create a universal sofl gap, g{oj) oc tu^, in the linear 
density of states. This soft gap is a manifestation of the sea-gull singularity— (see also 
Ref.-) predicted in the framework of the soft potential model for glasses. Recently it was 
shown that u;^ behavior of the DOS is indeed a universal feature in disordered systems 
for low frequency bosonic excitations which are not Goldstone modes^. 

The physical reason for this gap is very transparent. Due to anharmonicity there is 
always a blue shift of the soft oscillator frequency under the action of the static force 
/. For small cu this shift is proportional to |/|^^^ and it is anomalously large for small 
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forces. The density of states in the gap then can be estimated as follows 



g{uj) oc 



uj +00 

J ujidui j 5{to- af'^) df oc uj\ (1.2) 



As a result, the random static forces (together with anharmonicity) effectively "push out" 
oscillators from the low frequency range to higher frequencies, creating the universal soft 
gap, g{u) oc a;^. One can also see this gap in the context of mechanical stability of the 
system but of another kind. Due to local anharmonicity small frequencies cannot survive 
in the system in the presence of random static forces. In some sense they are not stable 
even though random forces cannot transform single-well potentials into double-well ones. 

The width of the w^-gap is of the order of the Boson peak frequency cot oc {SfY^^, 
where 6f is the width of the random force distribution P{f). The Boson peak frequency 
Ub is typically smaller™ than the characteristic frequency Uc determining the width of 
the vibrational instability region^ 



(1.3) 



The strong inequality Ub <C oJc occurs if Qoiujc) ^ Qoij^Jo) i-e. Qoiyj) is an increasing 
function of u (since Uc <^ ojq). As a result, giu) oc in the interval < a; < cuft. 
Frequencies higher than ujb are only weakly affected by the static forces. Therefore, the 
linear DOS, 5^(0;) oc uj (and the right wing of the Boson peak) which was created in the 
course of the vibrational instability is conserved in the interval ^6 ^ w < Uc- Since at 
small frequencies < i-jf, we have the universal cj^ gap in the density of states, the left 
wing of the Boson peak takes also the universal form, g{uj)/uj'^ oc cu^. As a result in the 
region of the vibrational instability, < < cjc? we have a universal behavior of the 
density of states g{ijj) {gioj) oc a;^ for < < tu;, and g{ijj) oc for a;^ < a; < oj^ with 
the Boson peak feature, independent of the initial DOS gQ{uj). At higher frequencies, 
above ujc-, we keep the initial DOS go{uj) almost undistorted. 

The Boson peak was the main topic of our previous papers^'^. In the present work we 
shall concentrate on the two-level systems, i.e. the level splittings due to the tunneling 
through the barriers separating the two minima of the two-well potentials. But this 
consideration is not independent of the Boson peak since we will see that the main 
parameters of TLS's will be strongly interrelated with the parameters of the Boson 
peak. Therefore, these two universal phenomena should be considered together. 

In the standard tunneling model the TLS's are often characterized by the so-called 
dimensionless tunneling strength C, Ref.-^ 

C = ^ (1.4) 

where P is the density of states of the TLS's, 7 the deformation potential, p the mass 
density of the glass and v is the average sound velocity. The experimental value of C 
for different glasses is small and varies in a narrow band between 10~^ and 10~^. In our 
theory such small numerical values for C will emerge in a natural way. 

Several authors^^i^Sii^i^^ proposed that the approximate universality and smallness 
of C in glasses may be a consequence of the interaction between the TLS's. Roughly 
speaking the idea was based on a mean-field approximation. The ith. TLS produces at 
a distance a deformation 

Si ~ (1-5) 
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n{E) c^-tl-c^l^. (1.9) 



where 7, is the deformation potential of the ith TLS. As the deformation is inversely 
proportional to rf the distribution function of the deformations in a glass is a Lorentzian 
with width 6e proportional to the total concentration N of the TLS's: 

Se ~ — -. (1.6) 

The energy Ei (the interlevel spacing) of each TLS changes under the deformation e as 

SEi = -fiS (1.7) 

and from Eq. (11.61) one sees that the energies E^ of the TLS's are distributed in the 
interval 6E: 

SEc^^. (1.8) 

For small energies, the density of states n{E) is independent of both the energy and the 
concentration of TLS's: 

N ^ fw'^ 

Je~~^ 

This is the result of a purely classical approach. In this approach the dimensionless 
parameter Cd 

C.==^.l (1.10) 

is of order unity rather than of the order of 10~^ — 10~^. This was the main difficulty 
of the theory outlined in Ref 's.-i^'^°'^^'^^. However, if one accounts for the quantum 
nature of tunneling, the situation is improved and the value of C is reduced strongly. To 
explain this on the qualitative level we remind that the energy £" of a TLS consists of two 
contributions, i. e. the classical asymmetry A and the quantum tunneling amplitude 
Aq: -E" = a/A^ + Aq. So far we have disregarded the latter. 
According to the standard tunneling model 

Aq = ^0 exp(— A), (1-11) 

where A is the tunneling parameter, distributed uniformly in the interval Amin < A < 
Amax- Usually Amin is taken to be about unity. Only TLS's with A ~ 1 can tunnel during 
typical experimental times. If Amax ^ 1, systems with A ~ Amax cannot tunnel and do 
not contribute to the observable properties. Therefore, the relative number of TLS's 
participating in the tunneling is proportional to the small number 1/Amax- As a result, 
we estimate the dimensionless parameter C in glasses as 

C ~ Cd/Amax ^ l/Amax- (1.12) 

If for example Amax — 10^ the dimensionless parameter C 10^'^. Thus the smallness 
of the dimensionless parameter C in our theory is related to typically large values of the 
tunneling parameter Amax (and to typically high barriers in the system). We will discuss 
this point in Section HVl 

In the same Section we will show that two important parameters, namely C for 
TLS's and the characteristic frequency Uc for HO's, marking the onset of the vibrational 
instability are closely interrelated: 
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Here W is an important characteristic energy in glasses^^. Typically it is of the order 
of a few Kelvin. It determines for example the position of the minimum in the reduced 
specific heat^iii, C{T)/T^ {W ~ 2Tinin) and some other properties of glasses above one 
KelviufSi. For vitreous silica ^ 4K. In particular it follows from this formula that 
the larger is the interaction I oc Uc between the oscillators the smaller is the TLS's 
parameter C: C oc 1//^. And as we will see in Section |IV] the smaller will be also the 
density of tunneling states, P oc 1//^. It naturally explains the very old puzzle in the 
physics of glasses, why the number of two-level systems is so small (one two-level system 
for a million of atoms). 

At a first glance this interesting result seems to be rather contradicting. The stronger 
the interaction / between the oscillators, the larger is the width of the vibrational in- 
stability region c^c oc / and therefore the higher is the number of double-well potentials 
created in the course of stabilization of the system due to anharmonicity. The explana- 
tion of this seeming paradox is that majority of the double- well potentials created due to 
vibrational instability have so high barriers V that they cannot participate in tunneling 
at all. As a result only a very small part of the double-well potentials contributes to the 
tunneling density of states, P. 

Since the experimental values of W and C are well known for many glasses one can 
estimate from Eq. (11.13^ the important characteristic energy hujc which gives the width 
of the vibrational instability region in glasses 

huj, ^ WC-^I^. (1.14) 

For example for a-Si02 W = 4 K and C = 3 ■ 10~^ giving hwc ~ 60 K. This falls into 
the Boson peak region {fkJb ~ 70K). As a result we see that indeed the Boson peak is 
placed in the vibrational instability range. 



II. VIBRATIONAL INSTABILITY 

To illustrate the idea of a vibrational instability, we consider a cluster containing a 
low-frequency harmonic oscillator (HO) with frequency ui surrounded by a large number, 
s — 1, of HO's with much higher frequencies Uj of the order of c^o ^ cui. Here uq is a 
order of magnitude estimate of the high frequencies. In glasses it usually is of the order 
of the Debye frequency. Let no be the total concentration of the HO's in the cluster and 
go{^^) the normalized initial density of states (DOS), i. e. the DOS of the HO neglecting 
their interaction. 

Including the interaction between the HO's, the total potential energy of the cluster 
is A; 1 1 

f/tot(a;i, X2, Xs) = ^ -^x^ ~ 2 ^ hjXiXj + ^ ^i^t^ > 0- (2-2) 

i * 

Here Xi are the generalized coordinates describing the vibrations of HO's, fcj > are the 
quasielastic constants of noninteracting oscillators and lij determines the bilinear inter- 
action between the oscillators. To stabilize the system we have added in this equation 
the anharmonic terms, Aix'^ (with Ai > 0). The interaction strength is given by^ 



(2.3) 
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where gij ~ ±1 accounts for the relative orientation of the HO's, r^j is the distance 
between HO's, p is the mass density of the glass and v is the sound velocity. 

The interaction between the HO's is due to the coupling between a single HO and 
the surrounding elastic medium (the glass). This HO-phonon coupling has the form^^ 

Hint = ^xe, (2.4) 

where A is the coupling constant and e is the strain. Introducing the masses of oscillators 
Mi we have for the bare frequencies (neglecting the bilinear interaction) as usual 



Vki/Mi. (2.5) 



These bare frequencies enter Eq. (12.11) for the initial DOS go{uj). We will assume the 
characteristic strength of the bilinear interaction / ^ Juq between the oscillators to be 
much smaller than the typical quasielastic constants, so that 

|/| < ^ ^ j^^^ (2.6) 

where M is the typical mass of the HO's. 

The equation of motion for the low-frequency oscillator is 

MiXi = -kixi + ^ hjXj - Aix\ (2.7) 
and for the high frequency ones 

For a slow motion one can set the acceleration term MjXj = in Eq. (12. 8p . For I <^ Muq, 
we have xj -C xi (see below) . Therefore we can neglect also the anharmonicity force term 
—AjX^ and the interaction terms IjiXi {i ^ 1) between the high frequency oscillators and 
get from Eq. 

Xj = (/,i/%)xi, J ^ 1. (2.9) 

According to Eq. (12. 6p we see that Xj ~ {I/Muq)xi -C xi. Inserting this value of xj into 
Eq. (12.71) we finally get a reduced equation of motion of the low-frequency oscillator 

MiXi = -{h - k)xi - A,xl = -'EL^!^ (2.10) 

dxi 

where 

^=y¥^Tr-2 (2-11) 

and the effective potential energy for the slow motion is 

f/eff(xi) = - t^)xl + ^A^xl. (2.12) 

The physical origin for this reduction to a one-mode approximation is the adiahatic ap- 
proximation where the high frequency modes adiabatically follow the slow low-frequency 
motion^^. As a result the interaction between the low and high frequency oscillators 
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renormalizes the quasielastic constant ki for the low frequency motion to the new effec- 
tive value 

k = ki- K. (2.13) 

For ki > K, the potential (12.121) is a one-well potential and the cluster of oscillators is 
stable, the equilibrium positions of all oscillators Xi = 0. For ki < k, the renormalized 
quasielastic constant k is negative and the cluster is unstable. The effective potential 
fl2.12p in this case is a symmetric double-well potential. This is what we call the vibra- 
tional instability. As a result of the instability the low-frequency oscillator is displaced 
to one of the two new minima 

xio = ±Vi^^h)lAi = ±^AWA'u (2.14) 

while the displacements of the high frequency ones are xjq = {Iij/kj)xio <^ xio (j 7^ 1) 
and are much smaller. The barrier height between the minima is 

As follows from Eq. fl2.12p the new lowest frequency of the system of s coupled 
oscillators is given by 

{ki~K)/Mi = k/M,, K<k,, 

^ \2(k- A;i)/Mi = 2|A;|/Mi, k > ki . ^' ' 

The first case k <k\ corresponds to a vibration in the minimum of a one-well potential 
fl2.12p while the second case k, > ki corresponds to a vibration in either of the two wells 
of a double- well potential (12.120 . It is remarkable that for weak interaction / <C MuJq 
the strength of the anharmonicity Ai does not enter the renormalized frequency (I2.16p . 

Using the Holtsmark method^^, we derived in our previous paper- the normalized 
distribution function of k 

IB ( 

P('^) = ^-?«exp -— (2.17) 



where 



Here Juq ~ /, and (1/00)^ ^ l/u^o is the u ^ moment of the normalized initial DOS, 
go{uj). This formula can serve as a definition of the important characteristic quantities 

ujc ^ //Mwo, and k^ = Mu^ = B^. (2.19) 

The physical meaning of these quantities is that the typical clusters with frequencies 
uji < Uc become unstable due to the interaction between the soft oscillator and the 
surrounding high frequency ones. Thus the characteristic frequency Uc indicates the 
onset of the mechanical instability region. We will see below that the creation of TLS's 
and the formation of the Boson peak occur in this region. 

However, in the present paper this particular form of the function p{k) is not suitable. 
The reason is the long-range power tail of this function, p(k,) oc 1/k^/^ for k, 3> kc. This 
tail leads to divergent integrals for large K-values when calculating averages of the type 
{k'^) for u > 1/2. As follows from Eqs. (12.110 and (12. 3p the long-range tail of the 



3V2Vm\cu/o ^ 
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distribution is related to close pairs with small distances between the low and high 

— 1/3 

frequency oscillators, Vij -C ng . However, usually the distance between the HO's in 
a glass can not be arbitrarily small and, therefore, the function p(/t) drops faster and 
approaches zero as k"^""'"^''/^ (for go{uj) oc cu", with n > 0). In the following we do 
not need the precise analytical form of this function. It will be sufficient to know that 
this function decays sufficiently rapidly for small and large k, with a characteristic scale 
K kc — I'^/ko- In Fig. [T] this function is shown for different interaction strengths J 
(no = M = uq = 1, see Ref.- for details). Knowing the function p{k.), we can calculate 



FIG. 1: Distribution function p(k), calculated as ensemble average by exact diagonalization of 
systems of = 2197 oscillators with go{uj) = 3u;^ and J = 0.07,0.10 and 0.15 (from left to 
right). 

the distribution function of the renormalized quasielastic constants, $(A;). Let F{ki) be 
a normalized distribution function of quasielastic constants fcj in Eq. fl2.2p . In the case 
of equal masses of the oscillators, Mj = M it is related to the normalized initial DOS 
go{uj) as follows 

F(k) = where k = Mu^ (2.20) 

and the normalized distribution function ^{k) is given as 

oo oo 

$(A;) = {5{k -ki + k))^^^^ = j dkiF{ki) j dnp{n)5{k - ki + k). (2.21) 



Integrating over the delta-function, it is convenient to present the expression for $(A;) 
for positive and negative k separately. We have from Eq. fl2.2ip 

oo 

$(A;) = j dkiF{ki)p{ki + |A;|), for A; < 0, (2.22) 
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and 

oo 

^{k) = j dKp{K)F{k + k), for A; > 0. (2.23) 



Since the distribution function p{k) is nonvanishing only for k < kc and rapidly drops 
to zero for k, > kc we conclude from Eqs. (12. 22^ and ( 12.23^ that the function $(/c) for 
\k\ <^ kc is approximately a constant 

oo 

<^{k) ^ $(0) = J dKp{K)F{K) ^ F{kc), for \k\ < fc,. (2.24) 



For negative k and \k\ > kc the function $(A;) rapidly drops. For positive k > kc, 
$(fc) ^ F{k). In Fig. [2] this function is shown for different interaction strength J 
(no = M = cuq = 1, see Ref.- for details). 



FIG. 2: Distribution function of the renormalized quasielastic constants, ^{k), calculated as 
ensemble average by exact diagonalization of a systems of = 2197 oscillators with go('^) = 
{F{k) = 3\/fc/2) and J = 0.07,0.10,0.15 (solid curves, from right to left). Dotted curve: 
Result of convolution, Eq. ([221]), for J = 0.1. 



III. INTERACTION BETWEEN LOW-FREQUENCY OSCILLATORS 

In the previous section we have considered the effect of the interaction between a 
low-frequency oscillator and the surrounding high frequency oscillators in a cluster. As 
a result of this interaction the quasielastic constant ki has been renormalized to a new 
effective value k = ki — n. Negative values of k indicate a vibrational instability of 
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the cluster. So far we have neglected the interaction between low-frequency oscillators 
belonging to different clusters. This interaction is much weaker and cannot produce 
a new instability. However, it causes internal random static forces acting on the low- 
frequency oscillators. As we have seen in the previous section, in the case of instability 
the interaction between low frequency and high frequency oscillators shifts the positions 
of their minima (static displacements). These shifts in turn act as forces if we take 
the interaction between unstable low- frequency HO's into account. As was shown in 
our previous paper- (see also Section |V]) these forces are responsible for the universal 
g{uj) oc u"^ (see also Ref.— ) dependence of the excess vibrational density of states for low 
frequencies. 

One can insert these forces into Eq. fl2.12p as a linear term —fxi where / is the 
internal random force created by the other unstable low-frequency oscillators. The 
effective potential energy then reads 

Ues{x) = -fx + hx^ + ^Ax\ (3.1) 

Here and henceforth the index 1 will be omitted. 

The distribution function of the random forces P{f) has been obtained in our previous 
paper- 

^(/) = ^7^. (3-2) 

where 5f is the width of the distribution. The Lorentzian form of the distribution 
is related to the fact that the forces between harmonic oscillators decay as [see 
Eq. ([53D]. 

One can estimate the width, 5/, of the distribution as follows. The static force 
exerted on the zth oscillator by the jth one is 

fi = lijXjQ. (3-3) 

Its characteristic value is given by displaced harmonic oscillators with frequencies of 
order of uj^- For these oscillators we have 

l[f ^ Jric, J ^ I/no ~ MuJcUJo/rio ^ MuJc/nogo{uJo). (3.4) 

Here Uc ~ nQgQ{uJc)uJc is the concentration of these unstable harmonic oscillators (double- 
well potentials) while no is the total concentration of HO's. Due to the normalization 
condition u}og{u!o) ~ 1. From Eq. fl2.14p it follows that the characteristic static displace- 
ment of these unstable oscillators is Xjo ~ = uJc\/M/A. As a result, one gets 
from Eq. (see Ref.-^) 



Sf^Jn^J'-L^MM.!^ (3.5) 

\ A \ A ^o(^o) 

As mentioned already in the beginning of this section, these internal random forces 
do not produce a new vibrational instability. However, they can transform some double- 
well potentials into single-well ones. For A; > the potential (13. ip is always one-well 
whereas for /c < the potential (13. ip is double- well for small forces |/| < where 

= 2\k\^^y3V3A. (3.6) 
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For I/I > (and A; < 0) (13.11) is a one- well potential. For \f\ = f^{k<0) the potential 
(13. ip is one- well with a bending point. 

It is interesting to compare the width of the distribution 6f with the characteristic 
value of They become equal for \k\ = k* = M {uj*)'^ where 

uj* ^ ujc [go{iUc)/goii^o)]^^^ < ^c- (3.7) 

The strong inequality u* ^ Uc occurs if both Uc ^ ujq and go{uJc) ^ go{uJo). We shall 
see in Section |V] that the frequency plays a role of the Boson peak frequency. 



IV. TWO-LEVEL SYSTEMS 

As follows from Eq. (13. ip negative values of k for / = correspond to symmetric 
double-well potentials. In a purely classical treatment the oscillator will vibrate in either 
of the wells. Taking quantum mechanics into account, there will be a finite probability of 
penetration through the barrier separating the two wells, i. e. there is a finite tunneling 
probability. This causes a splitting of the vibrational levels. We are interested in the 
lowest pair of levels. This constitutes a two- level system (TLS). These systems are 
ubiquitous in glasses and determine their low temperature properties^. 

Tunneling systems can be described effectively in terms of a tunnel splitting Aq and 
an asymmetry A. We will derive expressions for these quantities and their distributions. 
Neglecting the linear force term in Eq. (13.11) . the tunnel splitting is given in the WKB 
approximation as 

Xo XQ 

Ao^Wexp{-S/h) , S= j \p\dx = 2 j ^2M [U^^^^x) + V]dx. (4.1) 

—xo 

Here we approximated the dependence of the prefactor on the vibrational frequency by 
an order of magnitude estimate W (see also Ref.— ) 

W=^-{^\"' (4.2) 

2 V2M2 ) ^ ' 

that is of the order of the interlevel spacing in a purely quartic potential V{x) = Ax'^/4. 
From experiment, using the the soft-potential model^i^i^i^i^i^, one finds values W for 
different glasses of the order of a few Kelvin, e. g. for vitreous silica W ^ 4 K. The 
positions of the minima in the symmetric double- well potential Ues{x) (for / = 0) are 
denoted by ±xq = ±a/|A;|/A and V = k'^/AA is the barrier hight [see Eqs. (I2.14p and 

Evaluating the integral in Eq. f l4.ip we get 

^ ^ 2V2\kl^Hd^ 

and 



13 



The second quantity characterizing the TLS is the asymmetry of the two-well config- 
uration. For I/I <^ we have from Eq. (13. ip 



A = 2/xo = 2f^\k\/A. (4.5) 

We are interested in the two quantum states with the lowest energies. These states 
belong to both wells. TLS are often described by the interlevel distance E and the 
dimensionless tunneling parameter p: 



E=^Al + A2, p = {Ao/Ef . (4.6) 
The Jacobian of the transformation from the variables \k\ and f to E and p is 

d{E,p) \9j h? pVT^ ^ ' 

In the new variables the distribution function reads 

/oXVS ii^3/2tj/5/2 r -2/3 

F{E,p) = noPmmJl = [q) 25/'noP(0)$(0) (4.8) 

Here we have replaced P{f) and $(fc) by P{0) and $(0). This can be justified by 
estimating the relevant ranges of / and \k\ given by Eq. (14.111) . Using Eqs. (14.51) . (14.41) 
and (13.51) we express / in terms of Aq and A 

Taking rough estimates L = 10, E = 1 K, W = AK, huc = 100 K, uq/uJc = 3, and 
go{uj) oc cu^ one finds a typical value ~ 2 ■ 10"'^. To estimate \k\ we derive from 

Eq. (gaD 

|A;|/A;, = 4 (9/2)'/' L^/^ {W/hu,f . (4.10) 

This set of parameters gives the typical value of \k\/kc ~ 0.05. Therefore, in the range 
of parameters where the notion of TLS's is applicable, the characteristic values of / and 
\k\ satisfy the conditions 

/ < 6f, and \k\ < k^. (4.11) 

Our result can be compared with the standard tunneling modeli where the distribu- 
tion functions are 

P(A, Ao) = ^, F{E,p) = (4.12) 

Aq 2p^l-p 

with a constant density of tunneling states P. Comparing these distributions with 
Eq. (14. Sp . one gets 

P = (0 2'/'noP{0)m ^'^'^'^' L-'/^ (4.13) 

Both distributions, Eqs. (14. 8 p and (I4.12p . coincide regarding their dependencies on E 
and p, apart from the factor L~^/'^ describing a weak logarithmic dependence on E and 
p. The same factor is found in the soft potential model (see Refs.— ^^i). 
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To compare the tunneling strengths of the TLS's with experiment we study their 
interaction with strain, described by the deformation potential 7. According to Eq. (12 .4^ . 
the variation of asymmetry due to a strain e is 

= 2Kxq e = 2Ae^/\k\/A. (4.14) 
The deformation potential is defined as 

and from Eqs. (ITOj) . (531) and (g^D follows 

7 = Aa/V = / ^ • 4.16 

V A 27/6 y]^^ ^ 

In the standard tunneling model the TLS's are often characterized by the dimen- 
sionless tunneling strength C, given by Eq. (11.41) . For different glasses its value varies 
between 10~^ and 10~^. Using Eqs. (I4.13P and (14.161) one gets 



A^ W^/^^ 



C = 2V2noPiOMO) — (4.17) 



This value is independent of E and p as in the standard tunneling model. 

The different factors entering the expression for C can be estimated from our model 

as 

Aypv'' = J^I/no, I^Mio^Uo, ^0) ^ F{k,) = go{u;,)/2Mcu, (4.18) 
and from Eqs. (Q, and (g^) 

P(o).4-«i±:^^. (4.19) 

As a result we arrive at the important estimate 
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that is independent of the initial DOS of HO's go{^^) and in this sense is universal. It 
only depends on the characteristic energy W and the frequency Uc that is proportional to 
the interaction I. The larger the interaction between the original oscillators the smaller 
is constant C, C oc 

In the discussion of the Boson peak (Section V) we will see that Uc is two or three 
times larger than the Boson peak frequency Uf, (or u*) which slightly depends on initial 
DOS goiuJc) (see Eqs. (IS.lOp . (15.111) and Fig. [3]). Using values appropriate to Si02, 

= 4 K and huc = 100 K we get 

C ^ 10~^. (4.21) 



Thus the unified approach of this paper gives a value of the tunneling strength C in good 
agreement with experiment. Since values of C and W are well known from experiment 
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for many glasses^, Eq. (I4.20p can be used to estimate the important characteristic energy 
huc giving onset of the vibrational instabihty in glasses 

hu^ ^ WC'^I'^. (4.22) 

It was demonstrated in Ref.— that taking into account the experimental data for C and 
W this energy is indeed correlated with position of the Boson peak in glasses, cjf, ~ uj^. 
The two factors entering C, Eq. (11.41) . can be estimated separately as 

P^^fi^V^-^^ 4-— — ^^/^- (4-23) 

The first of these quantities is oc h}^^ and the second oc h?^^ and thereby C oc h. Since 
C ^ 1 is a dimensionless quantity it can be represented in the form C = h/S where 
5 ^ ^ is some classical action. From Eqs. f l4.20p . (14. 3 p and (14. 2 p we get 

on be o A 

i.e. S ~ SttSc- The classical action Sc corresponds to a typical double-well potential 
1 KB with \k\ = kc (and / = 0). 

Using the estimates uJc/ujq ^ 1/3, W = 4K and huc = 100 K we estimate the 
concentration of tunneling systems with energies in the range < E < W and tunneling 
parameter p ~ 1 as 

nxLS- W~no— ) ~ 3 . io-8no. (4.25) 

ujQ \nujcj 

The number of active TLS's is thus less than one for ten million of oscillators. This ex- 
plains why the concentration of observed TLS's in glasses is so small. Since according to 
Eq. (14.251) riTLs oc the number of TLS's decreases rapidly with increasing interaction 
strength /. 

It would be instructive to derive a dimensionless parameter Cc\ by a classical proce- 
dure (neglecting the tunneling probability Aq). We take the width of the force distribu- 
tion from Eq. (13. 5p and estimate the typical asymmetry Ac from Eq. (14. 5 p 

8f ^ I^'^^x^^ ^ Jn,^/kjA, A, ^ 5f^kJA ^ Jn, ■ {kJA). (4.26) 

With Uc the concentration of double-well potentials we get the classical estimate for 
their density of states 

P,.-.jL.-(iLV (4.27) 



Ac Jkc hujQ \hujc 

that is independent of h. 

For the deformation potential 7c we have from Eqs. (14.141) and (I4.15P the estimate 

..^A^, and 4«4«^(^)^ (4.28) 
also independent of h. Finally 

Cci = ^ ^ 1, (4.29) 

pv' 
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i.e. the dimensionless parameter Cd in this classical approach is of the order of unity 
which is a consequence of the interaction between the TLS'si^'^'^'^. We wish to 
emphasize that Cd (unlike C) does not determine any physical property of glasses. 

The reason for the difference between the two approaches (quantum and classical) 
is the following. In the classical approach we take all the double-well potentials into 
account. They have typically \k\ ~ kc- Their concentration, Uc, is unimportant since 
it is canceled in Eq. f l4.27p for Pc- In the quantum approach only the small portion 
of TLS's which are able to tunnel (they have \k\ << kc) contribute to the observable 
quantities. For all other TLS's the high barriers V and asymmetries A ~ Ac prevent 
tunneling (but they contribute to the internal random static force 6f). 

To further clarify this point let us consider TLS's with a quasielastic constant k in 
the interval Ak ~ k where \k\ <^ k* <^ kc [see Eq. (13. 7p ]. Their concentration n^, 
asymmetry A^ and the deformation potential 7^ are given by 

nk^no\kmO){rj6f), ^ r,./\k\/A, 7, = AvW4- (4-30) 

In this expression for we took into account, that for \k\ <^ k*, only the small fraction 
of all potentials, where -C 1, is of double- well type. Keeping in mind that $(0) ~ 

F{kc) we get the density of states 



and the parameter Ck 



Ak df^WJA 



Pkll noF{kc)J\kf/^ (\k 



3/2 



pv^ 6fVA \kcj 

Here we have used Eq. (14.261) for 5f and the estimate Uc — nokcF{kc) — noco'cfl'o('^c)- 

If we now fix \k\ by the condition that the exponent in Eq. (14.41) is of order of unity, 
i.e. \k\ ~ MW'^/h'^, we reproduce our quantum result, Eq. (I4.20p . Ck ~ (W/fuJcY- The 
classical result, C^, Eq. (I4.29P would be recovered for k ^ kc when ~ 1. We conclude 
that the physical reason for smallness of the parameter C for TLS's in glasses is the 
scarcity of those TLS's that are able to tunnel compared to their total number. 

Again we can compare our results with the standard tunneling model^. In this model 
the tunneling amplitude Aq = ^o6xp(— A) and the dimensionless parameter A is uni- 
formly distributed in the interval Amin < A < Amax- The lowest value Amm — 1- According 
to Eq. (14. ip \ = S/h, and therefore the maximal value Amax ~ Sc/ h. Taking into account 
Eq. fg:^ we get 

Amax ~ Sc/h ^ l/37rC. (4.33) 
Thus Amax is related to the small parameter C. For Si02, C = 3 ■ 10~^ and Amax ~ 350. 



V. THE BOSON PEAK 

In this section we relate the results obtained for the TLS's to the Boson peak prop- 
erties — see Ref.-. For this we calculate the vibrational density of states (DOS) g{uj). 
We start from the case / = (i.e. neglecting the interaction between the clusters). 
There are two types of harmonic vibrations. In the one- well case, A; > 0, according to 
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Eq. (12.241) the distribution function of k for /c -C /c^ is constant, $(fc) ^ $(0). Therefore, 
since according to Eq. (12.161) (top) k = Mu"^, the renormahzed DOS for u <^ Uc is 

gi{uj) = 2noM$(0)cu. (5.1) 

For harmonic vibrations in either well of a symmetric double- well potential [Eq. (13. ip . 

< and / = 0] \k\ = Muj'^/2 [see Eq. flTTBD . bottom] and ^{k) = $(0) for \k\ < k^ 
{uj < Uc) the DOS is 

giiiu) = noM$(0)o^. (5.2) 

It is half of the one-well contribution. The total DOS for / = and a; -C is the sum 
of the two contributions 

gtotiuj) = gi{uj) + gu{^) = 3noM$(0)o;. (5.3) 

It is a /mear function of uj independent of the form of the initial DOS go{uj). This linear 
behavior follows from the finite value of $(0). 

If the low- frequency HO's were isolated their density of states would be determined 
by Eq. (15. 3p . As we have shown in the Section IIIII there is, however, an interaction 
between these oscillators which we have to take into account. According to Eq. (13.30 
the low-frequency harmonic oscillators, displaced from their equilibrium positions (and 
forming the double- well potentials), create long-range random static forces / acting 
on other oscillators. In a purely harmonic case, these linear forces would not affect 
the frequencies. Anharmonicity, however, renormalizes the low frequency part of the 
spectrum, a manifestation of the so-called sea-gull singularity treated in detail in Ref.— 
(see also Ref.-). 

We begin with the case k > 0. It corresponds to one- well potentials. Consider 
an anharmonic oscillator under the action of a random static force /. The effective 
potential is given by Eq. (13. ip where ^/kjM is the oscillator frequency in the harmonic 
approximation for / = 0. The force / shifts the equilibrium position from x = to 
xo 7^ 0, given by 

Axl + kxo- f = 0, (5.4) 
where the oscillator has a new (harmonic) frequency 

Mul^ = k + 3Axl. (5.5) 

With $(/c) as the distribution function of k [see Eq. (I2.23p ] and P{f) as the distribution 
of random forces / [see Eq. (13. 2p ] the renormahzed DOS is given by 

oo oo 

gj{uj) = no [ Hk)dk [ dfP{f)6{uj-uJn,^). (5.6) 



Assuming a; -C and integrating Eq. (15.60 with $(A;) = $(0) we get the integral 

,M^2n,m^ f "k4Mm= (5.7) 







where according to Eqs. (15. 4p and (15. 5p 



1 / Af(j;2 — h 

f{k) = Axl + kxo = 3^ {2k + Moo'). (5.^ 
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Taking the Lorentzian distribution, Eq. (13. 2p . for P{f) and introducing a new variable 
t = a/1 — k/Muj'^ we finally get 

1 

= — noM$ — I — ) / — 5.9 



with 

to'' = VSA^/^SfY^^/VM. (5.10) 

The function gi{u) depends on a single parameter, u* characterizing, as we will see 
below, the position of the Boson peak Ub (cUb ~ uj*). The frequency u* is determined 
by the characteristic value of the random static force 6f acting on an HO with the 
characteristic frequency Uc. As a result, taking into account Eq. (13.51) . we get the 
estimate 

1/3 

uj''<^ujc. (5.11) 



UJ ^ LUr 



Again, as in Eq. (12.161) in lowest order the anharmonicity A does not enter this formula. 
This equation for u* coincides with Eq. (13. 7p obtained from the condition ~ 6f. 

According to Eq. (15.111) . for weak interactions / (cUc -C ujq), the frequency of the 
Boson peak u* <^ Uc only in the case when the initial DOS, go{^), is monotonically (and 
rapidly) decreasing (to zero) function of u. For example we can take goiuj) oc cu" with 
n > 0. Then for n = 2 and Uc = ujq/S we have from Eq. (15.111) u* ^ uJc/2. For the 
same n and smaller interaction, = ^o/^i we get uj* ~ ciJc/3. In the opposite case if 
the initial DOS drops to zero too slowly or remains nearly constant, go{uj) ~ const, we 
have from Eq. (15.111) that u* ~ ujc- In this case the Boson peak frequency uj* is of the 
same order as the characteristic frequency uj^. 

For small frequencies, uj <^ uj* , only small forces f Sf contribute to the integral in 
Eq. (15.71) . In this case the distribution function P{f) can be approximated by a constant 
value, -P(O), and we get from Eqs. (15. 7p and (15. 9p 



gi{uj) = 4no<l'(0)P(0) ^ = —noM^O)uJ (^— j oc uj\ (5.12) 



As a result, at low frequencies the renormalized excess DOS, gi{uj) oc uj^, Kei.—^. For 
sufficiently large frequencies, uj ^ uj* (but still smaller than uJc) the action of random 
static forces on the HO spectrum can be discarded. In this case the integral in Eq. ( 15. 9p is 
equal to {'k/6){uj*/uj)^. We recover the linear DOS, Eq. (15. ip . gi{uj) = 2noM^{0)uj oc uj. 

For k < the effective potential energy including a random static force / is given 
by Eq. (13.10 . The simple analysis in Section UTTl shows that for sufficiently small force, 
1/1 < where is given by Eq. (13. 6p . the potential (13.11) has two minima (double- well 
potential). For a large force |/| > the potential (13. ip has only one minimum (one- well 
potential) while for f = the potential is a one- well potential with a bending point. 

The calculation of the DOS for the lower minimum in the two-well case can be 
considered together with the one- well case (A; < 0, / > /*). The position xq of the 
minimum can be found from the equation 

Axl - \k\xo = f. (5.13) 

For / > and / < one should take, respectively, the positive and negative root of this 
cubic equation. In this minimum the oscillator has a harmonic frequency 

MujI„ = -\k\ + 3Axl. (5.14) 
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The density of states can be calculated from the Eq. (15.61) . For -C 

^„(o^) = 2no$(0)— = / d\k\-^M^M= (5.15) 



''iA J \/Mu^ + \k 
where according to Eqs. (I5.13P and (I5.14p 



Comparing f{k) with given by Eq. (13.61) one can see that the region of integration 
in Eq. (I5.15p . < \k\ < corresponds to the case of one minimum (one- well 

potential, f{k) > f^) and region < \k\ < corresponds to the case of two 

minima (double-well potential, f{k) < f^). Taking into account Eq. (13.21) for P{f) and 

introducing a new variable t = a/1 -l- \k\/Muj'^ we finally get from Eq. (I5.15P 



'3/2 

^jj(cu) = -noM<l>(0)— f— V / f-^ —. (5.17) 

As follows from this equation at small frequencies, u <^ u* , 9\\{.^) oc and for 
moderately high frequencies satisfying the inequality cu* <^ u; -C uJc, the integral in 
Eq. (I5.17P is equal to (7r/12)(ci;*/u;)^. Therefore in this case gii{uj) = noM$(0)ci; what 
coincides with Eq. (15. 2p . 

Combining results (15. 9p and (15.171) we get for the total DOS at T = 



/3/2 

g,^,(uj)=g,{uj)+gniuj) = -noMm-A^y [ ^ , , , Jl,,^ ^^2,2 - (5-18) 

vr u* Ku''/ J l + (ujuj)r(3-2rY 



We want to mention that Eq. (15.181) differs from the corresponding equation (22) of Ref .- 
not only by the prefactor and the upper limit but also by the power of (3 — 2t^) in 
the integrand. Thus we are correcting our error in Ref.™. Fortunately this does only 
marginally alter the plots of the quoted paper, where the analytical theory is compared 
to the results of simulation and experiment, since the plots of the present function 5'tot(i^) 
and the one given in Ref.- differ only slightly. 

For 00 <^ 00* the integral in Eq. (15.180 is equal to a/3/2 and we have 



(?tot(c^) = ^y|noM<l>(0)cu(-^j (xcu^ cu<cu^ (5.19) 

Taking into account Eq. (14. 2p one can present the density of states as function of energy 
E = hw (for E E* = huj*) in the following way 

n{E) = gUu^)/h = —nomPiO) ^ j • (5-20) 

This result for n{E) can be compared with Eq. (14.80 giving the density of states for 
TLS's. It is clear from this comparison that for E ^ W the density of states of HQ's 
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is much bigger than the density of states for TLS's. For large frequencies uj* uj <^ uj^ 
the integral in Eq. (15.181) is equal to (vr/4) ■ {uj* /uY and we have 

5'tot(t^) = 3noM$(0)u;, cj* < < cj^, (5.21) 

which coincides with Eq. (15.31) . 

Since at low frequencies uj <^uj* , the total DOS gtoti.'^) oc a;^ and at high frequencies 
cu* <^ -C ujc, Qtotil^) oc cij we have a peak in the reduced density of states gtot{(jj) / u)'^ 
at u; ~ u;*, the Boson peak. In the Fig. [3] we plot the function gtotij^) / oj'^ ■ We see from 
this figure that uo^ oj*. This figure is valid only for the case uj* ^ uoc- If oo* — ooc 



FIG. 3: The Boson peak: the reduced density of states (7tot(w)/a;^ given by Eq. (jS.lSp . 

{goiuj) ~ const) then uj^ ~ ojc and at u; ujc, according to our previous results (see 
Ref.-, Eq. (27)), gtot{oj) ~ goioj) ~ const. In this case the right wing of the Boson peak 
is determined by the initial density of states gQ{uj) and gtotiyo) / oj"^ oc l/oj"^ (instead of 
gtotjuP' oc 1/u; in the previous case). 

The DOS for the higher minimum in the double-well potential (13.11) is different. 
Though the thermal occupation number of this minimum is smaller than the one of the 
lower minimum it can contribute to the total DOS at finite temperatures. Starting from 
Eq. (15.61) the position of the higher minimum can be obtained as the smallest negative 
root or largest positive root of Eq. (15.131) for / > and / < 0, respectively. The resulting 
DOS for a; <C is 

with 



f{k) = Axl - \k\xo = y ^''3^'^' (2|fc| - Mco'). (5.23) 

Inserting the Lorentzian distribution for P{f), Eq. (13. 2p . and introducing a new variable 
t = v^l + \k\/Mu^ we get for the DOS 

oo 

g^^^[uj) = -noM^O)—(—y [ f-^ —. (5.24) 



3/2 
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At low frequencies, uj -C cu*, the integral in Eq. (15.241) is equal to (vr/S ■ 2^/^)(c<j*/u;) 
and ^ 

giiiiyj) = 2^/^noM$(0)cj^ (—) oc < cu^ (5.25) 



We see that at low frequencies the dependence of gm^uj) (Eq. fl5.24p ) differs from u;^ 
dependence, the DOS in the higher minimum is proportional to uj^ in accordance with 
Ref.^. From Eqs. (I5TT9|) and ([ESSD it follows that giu{u) / gu{uj) ^ u*/uj > 1 and 
therefore, for u u* and equal population of both minima the DOS in the higher 
minimum is larger than the one in the lower minimum. However, as we will see below, 
including the thermal population factors for the two minima reverses this. 

For high frequencies, a; 3> a;*, the integral in Eq. fl5.24p is (7r/12)(u;*/a;)^ and 

gui{uj) = noM^{0)uj, u -> (5.26) 

which coincides with Eq. (15.21) . 

So far we disregarded the thermal population factor. Taking it into account we get a 
temperature weighted DOS 

oo 

gM-2n.m^ I d\k\4M^ L__ (5.27) 

^ ^ ' ' ^ ' VSA J ' ^ ^Muj^ + \k\l + exp{A/T) ^ ' 

where A is the energy difference between the minima. 

In the low-frequency case, uj <C a;*, the integral in Eq.. (15.271) is a constant plus some 
cu-dependent correction. To estimate the constant, we set u; = in the integral. The 
higher minimum then turns into a bending point and 

oo 

/ ,M2o;3 [ d\k\ Pin) 
~giii{LO,T) = 2no$(0)— = / (5.28) 



2>A J ^l + exp(A/r)' 
^ 

where A = 'ik'^/AA is the energy distance between positions of the minimum and the 
bending point in the potential Uce{x) for f = fl [see Eq. (13.61) ]. 

The result of the integration in Eq. (15.281) depends on which of the two functions 
P{fk) [1 + exp(A/T)]~ decays faster with \k\. The function P{fk) decays with a 
characteristic scale 

\k\ =kf = 3Ai/3(5/)2/3/22/3 (5.29) 

(for = Sf). The function [1 + exp(A/T)]~^ decays with a characteristic scale \k\ = 

kr = 1\/ AT I'i. Both scales become equal at the temperature T = T*, where T* is given 
by 

^ - 8:2V^^4V^ - 8:2V^^^'^ ^ -T28:2V^^ \W ) ■ ^^-^^^ 

Estimates show that T* is a rather large. For example for fiuj* = 40 K and = 4 K we 
have T* ^ 740 K. 

Therefore, the low-temperatures case is more realistic. For T ^ T*, one has kx kj 
so that P{fl) = P(0) and therefore from Eq. (Km 

oo 

9ni(-, T) = ||„„*(0)P(0)^rV* I « (5.31) 
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The last integral in this equation is equal 1.07 ~ 1. Therefore the equation can be 
rewritten in the form 

^,„(,,T)«^^„„t(0)M(-) -(-) . (5,32) 

Now let us compare Qtotil^) and ^iii(a;, T) for oj <^ . Taking into account Eqs. (15.191) 
and f l5.32p we have 



Thus for T <C {hj/W)'^ we get 5'tot(t^) ^ gm{^)- In the opposite case the contribution 
of the higher minimum to the DOS dominates. 



VI. RESONANT SCATTERING OF PHONONS BY HO's 

Taking Eq. (12.41) for the coupling of the quasilocal oscillators to the phonons we get 

7-1 '^^'^ 

^-ho = 2m;^^(") ^^-'^ 

where /res,HO is the mean-free path of phonons due to resonant scattering on quasilocal 
HO's with a density of states gi^uj). For low frequencies, below the Boson peak frequency, 
uj ^ u*, we have from Eq. (15.201) (see also Ref.— ) 

where C ~ 10^'^ ^ 10"^ is the TLS's dimensionless parameter given by Eqs. ( 11. 4p and 
(I4.17P (see also the estimate, Eq. (14.201) ). Its value is well known from the low temper- 
ature properties of glasses. For high frequencies, above the Boson peak, in the interval 
Col* -C -C cUc, we have from Eq. (I5.2ip . gioj) = 3noM^{0)u. As a result 

Cho = 5T^*(0)^«- (6.3) 

Let us compare the last quantity (proportional to uj) with the inverse wave length of 
the phonons = uj/2ttv. We have the ratio 

^ 37r^^<|.(0). (6.4) 



/res,HO P^^ 

Using the estimates 

K^no/pv^ ^ / ^ Muo,uq, $(0) ^ F{K) = g^{uj,) /2Muo,, (6.5) 



and ojq ^ l/(7o(u;o) we have 



A ^3^,^^3^, /^y^ 



^res, HO 2 goiijQ) 2 \uj, 
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the last estimate follows from Eq. (15. lip . Thus this ratio is a constant in the interval 
uj* < uj < ujc and depends only on the characteristic frequency Uc oc I and the behavior 
of the initial DOS go{uj). From the last equation it follows that ratio A//res,HO depends 
on the cube of the ratio of two important frequencies, the Boson peak frequency uJt ~ uj* 
and the characteristic frequency Uc- Both of them can be measured on experiment (see 
Eq. fl4:22|) and Ref.i^). 

For the weak interaction / which we consider in the paper Uc ^ u^o and if for uj ^ 
initial DOS go{uj) also goes to zero sufficiently rapidly, then go{uJc) ^ fl'o('^o) and A -C 
^res,HO; i-G. resouaut phonon scattering is also weak. However, due to big numerical 
coefficient in Eq. (16. 6p in some realistic cases we can have a strong phonon scattering. 
For example if the initial DOS go{uj) oc u"^ and Uc ~ ujo/3, then 

A//re.,HO ~ 7rV6 ^ 1.64. (6.7) 

In this case criterion of loffe-Regel for the phonons (/res.HO < A) is approximately satisfied 
and we have strong phonon scattering above the Boson peak frequency in the interval 
UJ* u <ti ujc (which is not too big since in this case uJc ~ 2uj*). 

Another interesting case is a flat initial DOS, go{uj) ~ const, then go{uJc) — goiuJo) 
and as follows from Eq. (15.111) uj* ~ Uc and interval [ci;*,to'c] shrinks to one point (Boson 
peak frequency, ujh ~ uJc) and at the Boson peak frequency we have 

= - 15. (6.8) 

In this case the criterion of loffe-Regel is again satisfled and we have at the Boson peak 
frequency uj* ~ uJc the regime of very strong resonant scattering of phonons on quasilocal 
harmonic oscillators, independent of the strength of interaction / (and uJc). However, 
it is necessary to stress that in this case the strong scattering takes place only in the 
vicinity of the Boson peak frequency uJc- 

At higher frequencies the initial DOS go^uj) ~ const and according to Eq. (16.10 
the resonance phonon mean free path is also constant and independent of frequency. 
However the phonon wave length A oc 1/u; decreases with frequency. Therefore the 
regime of the weak phonon scattering will recover again at higher frequencies uj ^ uJc- 
Similar behavior was observed in Ref.^^ for resonant scattering of phonons on librational 
(quasilocal) modes in crystals. We can extend the regime of strong scattering to well 
above the Boson peak frequency (up to Debye frequency ujq) only when the initial DOS 
is a linear function of the frequency, go{uj) oc uj and interaction / is not too small, 
UJc > cfo/15 (see Eq. (16. 6p ). 

As we already mentioned, we will have a weak resonant scattering of phonons on 
quasilocal oscillators only when the interaction / is sufficiently weak and the initial 
density of states go{uj) decreases to zero sufficiently fast with uj, so that go{uJc) <^ 
(2/37r^) (7o(^o)- In such a case the mean free path of the phonons /res.HO "will be much 
larger than their wave length A in the whole frequency range. In this case phonons are 
well defined quasiparticles everywhere. 

We give here also the relaxation time r of a HO with frequency uj due to the interaction 
with phonons. From Eq. (12.40 we get 

1 _ A=^ uj^ _ JW 

Estimating 

A ~ So/a, pv'^ ^ ^o/a^ Mv"^ ~ £o, fiv/a ^ /zt^o, (6.10) 



(6.9) 
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where — 10 eV is of the order of atomic energy, a ~ 1 A is of the order of interatomic 
distance, and uq is of the order of Debye frequency we get 

l/ujT ^ lu/Attluo <^1. (6.11) 

Therefore HO's with u <^ uq are well defined objects. 



VII. DISCUSSION 

In our previous papers^i^ we proposed a mechanism of the Boson peak formation. 
The essence of the mechanism can be formulated as follows. A vibrational instability 
of the weakly interacting QLV's (stabilized by the anharmonicity) is responsible for 
the Boson peak in glasses and other disordered systems. The instability occurs below 
some frequency Uc proportional to the strength of the interaction I between low and 
high frequency oscillators. Whereas anharmonicity is essential in creating the atomic 
structures supporting the Boson peak, the vibrations forming the peak in the inelastic 
scattering intensity or the reduced density of states are essentially harmonic. 

The present paper extends these ideas. We show that such seemingly unrelated phe- 
nomena in glasses (typical for the glassy state and usually treated by separate unrelated 
models) as the formation of the two-level systems and the Boson peak in the reduced 
density of low-frequency vibrational states g{(jj)/uj'^ can be explained by the same phys- 
ical mechanism, namely the vibrational instability of weakly interacting soft harmonic 
vibrations. These can be seen as localized vibrations with a bilinear interaction with 
the extended modes, the sound waves. The resulting exact harmonic eigenmodes are 
quasilocalized vibrations that have been observed in numerous computer simulations 
and have been discussed extensively — see Ref.- and the references therein. 

The instability, which as in all solids is controlled by the anharmonicity, creates a 
new stable universal spectrum of harmonic vibrations with the Boson peak feature as 
well as double-well potentials with a wide distribution of the barriers heights that is 
determined by the strength of the interaction I between the oscillators. Depending 
on the barrier height (and temperature) these will lead to tunneling and relaxational 
transitions. To check for the consistency of our theory we calculated the dimensionless 
parameter C = P'j'^/pv'^ ^ 10~^ for the two-level systems in glasses which is observed in 
experiments^. The smallness of this parameter is a longstanding puzzle. In our theory 
it follows naturally. The physical reason for small value of the parameter C is that only 
a small fraction of all created TLS's can actually tunnel in realistic timescales. 

We show that the larger is the interaction / between the original harmonic vibrations 
the smaller is parameter C. It reminds partly the ideas of Ref.= though we do not have 
here the frustrated strong interactions. We prove that for our simple model weakly in- 
teracting oscillators C = (W/hjc)^ oc Here W is an important characteristic energy 
in glasses of the order of a few Kelvin. The value of C is independent of the assumed 
initial DOS of HO's go{^^) and in this sense it is universal. Varying the characteristic 
energy W and interaction I (i.e. the characteristic energy huc) for different glasses C 
lies in the interval from 10~^ to 10""^. However, we want to stress that we are not free in 
the choice of these two parameters. The energy W is well known from experiments on 
specific heat^iii, thermal conductivity^^ and heat release^ in glasses. As for the charac- 
teristic frequency Uc it should be of the order or larger than the Boson peak frequency 

UJb. 

In the unified approach developed in the present paper the densities of tunneling 
states and of excess vibrational states at the Boson peak frequency are interrelated. 
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Since the experimental values of C and W are well known for many glasses we can 
use this formula to get the important energy fkUc = WC~^^^ giving us the onset of the 
vibrational instability region. For vitreous silica hwc ~ 60 K falling perfectly into the 
Boson peak range. The same holds for many other glasses^. It indicates that the Boson 
peak is indeed placed inside the vibrational instability range. 

It is instructive to compare the results of the present paper and of Refs.->^ with our 
earlier paper— where we also discuss the possible origin of the Boson peak. In this 
paper we consider low-energy Raman scattering in glasses. As in the present paper 
and in Refs.^^'i we assume that the scattering and the energy transfer are due to the 
interaction of the light with the soft potentials in glasses. The density of states of the 
quasilocalized HO's, according to Ref.— , is proportional to u'^ for low frequencies and 
to u for high frequencies. This behavior qualitatively resembles the one obtained in 
Ref.^. However, the phenomenon of vibrational instability was disregarded in-- and the 
discussion of the Boson peak was necessarily somewhat qualitative. Considering the 
vibrational instability, puts the theory on a more quantitative level and, for instance, 
allows the determination of the shape of the Boson peak. 

Let us now compare the results of our paper with previously published important 
class of models of the Boson peak— i^i^i^. Experiment has shown that the Boson peak 
is formed by largely harmonic vibrations. Therefore, in all these models the authors 
have considered an Hamiltonian of the form 



with randomly distributed quasielastic constants kij. Since this potential energy is purely 
harmonic we call all such models of the Boson peak harmonic random matrix (HRM) 
models. The main difference between the quoted four HRM models is in the different 
distributions of the quasielastic constants kij. 

If all quasielastic constants are positive, k^j > 0, then the corresponding dynami- 
cal matrix (Hessian) is positive-definite and therefore all the eigenvalues are obviously 
positive as well, utf > (z = 1, 2, ...n) excluding those zeroes which come from the trans- 
lational and rotational invariances. In such a case the system is mechanically stable. 
As was shown in^iiS^i^i the system remains to be stable even when some (rather small) 
fraction of the quasielastic constants kij in Eq. (17. ip is negative (and small enough). In- 
creasing the fraction of negative kij (or their absolute values) authors have approached 
the mechanical stability threshold (when the first negative u"^ has appears in the spec- 
trum) . 

The reduced density of states g{uj)/uj'^ for Hamiltonian (17.11) usually has a maximum 
at some frequency cUmax for typical values of kij > 0. Changing the parameters of the 
distribution function P{kij) one can shift this maximum to higher or to low frequen- 
cies. The last case was the main goal of the papei^^^^. The biggest red shift of the 
maximum has been achieved approaching the mechanical stability threshold. In the first 
two papers^i^^ the original maximum was due to Van-Hove singularity of the crystalline 
DOS. Atoms in these papers were placed on a perfect cubic lattice. As a result the 
Boson peak has been ascribed to the lowest Van Hove singularity shifted due to disor- 
der. In another two papers^^i^i atoms were distributed randomly in 3-d space (so-called 
Euclidean Random Matrix models) and quasielastic constants kij depend exponentially 
on the interatomic distances |rj — rj|. Therefore the distribution function of quasielastic 
constants in these two cases P{k) oc 1/k and has a singularity for k ^ 0. Thereby the 
portion of small k was increased compared to Gaussian and box distributions used in 
the papers'^^'^. Obviously the red shift of the original peak in g{uj)/uj^ was much more 
pronounced in these latter two case^'^. 




(7.1) 
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Our approach, which continues from our previous papers^i^ differs essentially in two 
ways. First, we postulate that the excess in vibrational modes originates from quasi- 
localized vibrations. Their existence has been shown in numerous simulations of different 
types of materials. Such modes can be described as local modes (cores) which weakly in- 
teract bilinearly with the extended modes (sound waves) and thus with each other. The 
exact harmonic eigenvectors are of course extended as in HRM models. Secondly, we do 
not invoke special distributions for the elements of the dynamical matrix to avoid un- 
stable vibrations but, on the contrary, show that the generic instability, when controlled 
by the anharmonicity which is present in all real systems, automatically gives both the 
TLS and the Boson peak (with a universal shape) without any further assumptions. 

The essence of the mechanism can be formulated as follows. The randomly distributed 
weakly interacting QLV become unstable at low frequencies in harmonic approximation. 
This is the equivalent to the instability in the general HRM models. Instead of assump- 
tions on distribution functions of interactions lij we use the always present anharmonicity 
as the stabilizing factor. The previous vibrational instability of the weakly interacting 
QLV's thus becomes responsible for the Boson peak and TLS's in glasses and other dis- 
ordered systems. Whereas anharmonicity is essential in creating the atomic structures 
supporting the Boson peak, the vibrations forming the peak in the inelastic scattering 
intensity or the reduced density of states are essentially harmonic. Comparing with ran- 
dom matrix models our Eq. (12. 2p without the anharmonicity term would correspond to 
the case of the unstable random matrix whereas the result of Eq. fl2.16p correspond to 
the stabilized case. And, importantly, the anharmonicity strength A (thanks to mirror 
transformation) does not enter in the expression for the renormalized frequency. So, 
the anharmonicity reconstructs the spectrum but the final result is independent of the 
strength of anharmonicity. The advantage of our approach is that the stabilization is 
not a result of an additional assumption but is a benefit of the vibrational instability + 
anharmonicity which is haunting the alternative approach. 

Summarizing briefly we can say that the Boson peak in papers2ii2^i2^i2i was obtained 
by a purely harmonic ansatz inside the mechanically stable region of their Hamiltonians. 
The position and the form of the peak depend strongly on the distribution function 
of quasielastic constants P{kij). In our approach the Boson peak is built inside the 
mechanically unstable region of the harmonic potential parameters. Therefore, the role of 
anharmonicity as stabilizing factor is crucial. But the form of the Boson peak appears to 
be universal and independent of the initial assumptions about interaction or the initial 
distribution function go{uj) or of the anharmonicity strengths Ai. Relating TLS's and 
Boson peak parameters in our theory we are able to show that the instability crossover 
frequency Uc given by Eq. (11.141) lies in the Boson peak region. Due to weakness of 
the interaction I the universal reconstruction of the spectrum in our theory takes place 
in the low frequency range only, leaving high frequency range nearly unchanged. This 
is different in the HRM models^i^i^i^ where the whole spectral range is completely 
reconstructed in the course of the change of the distribution function P{kij). 

As mentioned before in our theory the estimate for the dimensionless tunneling 
strength C emerges in a naturally (see Eq. (I4.20p ) . Varying the characteristic energy W 
and the interaction strength / for different glasses C falls into the interval between 10~^ 
and 10~^. We want to mention that other explanations of the smallness of the constant 
C have been proposed. Burin and Kagan^^ predict a number of universal properties 
of amorphous solids, including the small values of the tunneling strength C, due to a 
special form of interaction of such defect centers with internal degrees of freedom like the 
TLS's. It is important that the interaction between such centers falls off with distance 
r as In this case effects of correlations between many TLS's {dipole gap effects) 

might be important at sufficiently low temperatures. 
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In our paper we have neglected all such many particle correlation effects between 
tunneling TLS's which could further reduce (or stabilize) the value of C. One should, 
however, keep in mind that the law of interaction is valid as far as the effects 
of retardation play no role. In dielectric glasses these effects are determined by the 
sound velocity. At sufficiently large distances they could play a role that would result in 
the variation of the interaction law. We believe that the role of the retardation effects 
deserves a special investigation. 

Another important difference between our and Burin and Kagan approaches is that 
the parameter C in our theory is a quantum mechanical quantity, C oc h (see Eq. f l4.24p ). 
It disappears in the classical limit — > 0. In other words the smallness of C is directly 
related to the smallness of the quantum mechanical probability for particles tunneling 
through high barriers in a glass (see also Eq. fl4.33p ). Differently, the dipole gap effect 
is purely classical since it is based on the classical dipole-dipole interaction between 
TLS's. Therefore, it would be very interesting to elucidate which of two mechanisms 
(or both) dominates and is responsible for the small value of parameter C in glasses. In 
particular it would be very interesting to calculate the relaxation times for many particle 
correlations effects to build the dipole gap. 

An analysis of the low temperature properties of glasses along the lines of the soft 
potential model based on a numeric simulation of a Lennard- Jones glass was presented 
by Heuer and Silbey^S. They numerically searched for the energy minima of the glass and 
constructed double-well potentials for close minima. By extrapolation to small values 
they were able to extract the distribution functions for the soft potential parameters. 
These potentials correspond to our potentials after inclusion of the interaction between 
the HO's. The present theory is in agreement with their simulation results. 

The theory presented in this paper deals with the effects of soft modes produced by 
disorder which can be expected to have a broad frequency distribution. In the literature 
the term Boson peak is rather loosely defined. Often it is used for any low frequency 
maximum in the reduced DOS. In particular in plastic crystals, see eg. Ref.— , soft HO 
are present even before disorder. Consequently disorder only broadens their sharp DOS. 
Depending on the strength of this broadening our theory will apply more or less to such 
cases. The same applies to TLS which also can be present before disorder. 

In summary, we have shown that the same physical mechanism is fundamental for 
such seemingly different phenomena as formation of the two-level systems in glasses and 
the Boson peak in the reduced density of low-frequency vibrational states g{uj)/uj'^. In 
this way two of the most fundamental properties of glasses are interconnected. 
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